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Comparison of AGN and Nuclear Starburst Activity in Seyfert 1 and 2 Galaxies 
over a Wide Luminosity Range Based on Near-infrared 2-4 fim Spectroscopy 



O 

(N 

-t— > 
O 

O 

in 



< 
O 

6 



Nagisa Ol 1 ' 2 ' 3 
nagisa. oi @ nao. ac.jp 
Masatoshi IMANISHI 2 ' 1 ' 3 & Keisuke Imase 1 ' 2 
1 Department of Astronomical Science, The Graduate University for Advanced Studies (Sokendai), Mitaka, Tokyo 181-8588, Japan 
2 National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 
3 Present Address; Subaru Telescope, 650, North A'ohoku Place, Hilo, HI 96720, USA 

(Received ; accepted ) 
Abstract 

We present near-infrared K- (1.9-2.5/im) and L- (2.8-4.2{im) band spectroscopy of 22 Seyfert nuclei. We use 
two methods to investigate the presence of nuclear starbursts: (1) the Polycyclic Aromatic Hydrocarbon (PAH) 
emission feature at X res t = 3.3 /im in the rest frame of L-band spectrum (a starburst indicator) and (2) the CO 
absorption feature at X res t = 2.3-2.4 /im in the rest frame of the K-band spectrum, originating in the CO molecule. 
We clearly detected the 3.3 /im PAH emission features in five objects and the CO absorption features in 17 objects. 
Seyfert 2 galaxies tend to show bluer K — L colors compared with Seyfert 1 galaxies. We interpret the discrepancy 
as resulting from relative strength of stellar emission because AGN emission is affected by dust extinction. The 
3.3 /jm PAH emission luminosity (L^^pah) distributions for the Seyfert Is and Seyfert 2s are very similar when 
normalized to the AGN power. Star-formation rates estimated from L^^pah could be large enough to inflate the 
dusty torus by supernova explosion. We find that L^^pah positively correlates with N-band luminosity with small 
aperture over a wide AGN luminosity range, and is independent of physical area we probed. The results suggest that 
nuclear region has a concentration of star formation and the star formation would control AGN activity. 

Key words: active — galaxies:nuclei — galaxies:Seyfert — galaxies:starburst — infrared:galaxies 
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Introduction 



A galaxy that exhibits bright emission in the nuclear re- 
gion is termed an Active Galactic Nucleus (AGN). The ori- 
gin of the radiative energy is putative release of gravitational 
energy through accretion of the interstellar medium from the 
host galaxy onto a central super-massive black hole (SMBH) 
with a mass of >10 6 M Q . Starbursts are one of the candidates 
for transporting material onto the center (Wada & Norman 
2002), so studying the correlation between AGNs and star- 
bursts is crucial to understanding the origin of AGN activity. 
Seyfert galaxies are the second most numerous AGN class in 
the local universe 1 . They are classified into two types: type 1 
Seyfert galaxies show broad and narrow optical emission lines, 
and type 2 Seyfert galaxies show narrow emission lines only. 
The difference is explained by the AGN unification paradigm 
( Antonucci 1993) in which molecular gas and dust with a torus- 
shaped structure (the so-called dusty torus) surrounds the cen- 
tral SMBH. The intrinsic AGN properties of the two types of 
Seyfert galaxies are assumed to be the identical. The dusty 
torus plays various roles in observed spectra: (1) obscures the 
AGN for a Seyfert 2 galaxy and (2) absorbs the energetic ra- 
diation from the AGN and re-radiates it as infrared dust emis- 
sion for both types of Seyfert galaxies. Because the torus is 
composed of significant amounts of molecular gas and dust, 
it is a natural site for star formation. In particular, the outer 



portion of the torus is thought to be an ideal site for star for- 
mation because star formation can be suppressed in the inner 
portion by tidal forces from the SMBH. In fact, it was argued 
that such starburst activity could occur in the dusty torus in the 
nuclear region ( < lOOpc) of a Seyfert galaxy (Wada & Norman 
2002; Imanishi 2003). We call such starbursts in Seyfert nuclei 
"nuclear starbursts". Wada & Norman (2002) showed that the 
nuclear starbursts in the dusty torus can produce an inflated tur- 
bulent torus around the central SMBH. However, with limited 
spatial resolution, it is difficult to detect the emission from the 
nuclear starbursts as they are diluted by strong AGN emission. 
Here, we discern these emission and estimate the magnitudes 
of starburst activities using infrared K-band (1.9-2.5/im) and 
L-band (2.8-4.2/im) spectra. 

Polycyclic Aromatic Hydrocarbon (PAH) emission at 
around 3.3 /im in the L-band is a powerful tool for distinguish- 
ing starburst emission from AGN emission. PAH molecules are 
widely distributed throughout interstellar space (Tanaka et al. 
1996). The PAH molecules are not destroyed and are excited 
by non-ionizing UV photons from stars and then emit the line 
at 3.3 /im in the Photo-Dissociation Region (PDR) for starburst 
(Sellgren 1981), whereas the PAH molecules are destroyed by 
X-rays (Voit 1992) from AGN. Observationally, a normal star- 
burst galaxy, which consists of HII regions, PDRs, and molec- 
ular gas (e.g. M82 and NGC253), shows the 3.3 fim PAH 
emission (Tokunaga et al. 1991; Imanishi & Dudley 2000), 
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whereas a pure AGN shows only PAH-free continuum emis- 
sion from larger-sized AGN-heated hot dust grains (Moorwood 
1986; Imanishi & Dudley 2000) because an AGN emits strong 
X-rays in addition to UV Additionally, because the 3.3 /zm 
PAH emission feature is very strong, it is detectable even if the 
starbursts are weak. 

Nuclear starbursts in Seyfert galaxies can also be inves- 
tigated through infrared K-band spectra. The CO(Av=2-0) 
absorption feature at X rest =2. 3-2.4 /im is caused by the CO 
molecule in the photosphere of cool stellar population (i.e. red 
giants and supergiants) (Frogel et al. 1978; Doyon et al. 1994). 
Empirically, starburst galaxies exhibit a substantially deeper 
CO( Av=2-0) feature than do quiescent spirals galaxies (Frogel 
et al. 1978; Ridgway et al. 1994), and AGNs do not exhibit this 
absorption feature. 

Furthermore dust extinction is much lower in the K- and L- 
bands than at shorter wavelengths, Ak ~ 0. 062^4 y and Al ~ 
O.QSlAy, respectively (Nishiyama et al. 2008; Nishiyama et 
al. 2009), so that the quantitative uncertainty of the dust extinc- 
tion correction is significantly reduced compared with shorter 
wavelengths (e.g. UV, optical). Thus, K- and L-band spec- 
tra enable us to estimate the relative contributions of AGN 
and starbursts emission to an observed spectrum in a quanti- 
tatively reliable manner (Oliva et al. 1999; Ivanov et al. 2000; 
Rodriguez- Ardila & Viegas 2003; Imanishi 2003; Imanishi & 
Alonso-Herrero 2004; Imanishi & Wada 2004; Davies et al. 
2007). 

The above authors studied mainly nearby Seyfert galax- 
ies with high infrared luminosities (L/# >10 44 ergs s~ 4 ), and 
found a correlation between the AGN activity and nuclear 
starburst activity (Imanishi 2003; Imanishi & Alonso-Herrero 
2004; Imanishi & Wada 2004). However, it is unclear whether 
the correlation holds over a wide range of AGN luminosity. 
Haas et al. (2005) argued that the height of the torus is smaller 
in lower luminosity Seyfert galaxies compared with higher lu- 
minosity Seyfert galaxies, suggesting that nuclear starbursts 
are weaker in low-luminosity AGNs. In contrast, Ballantyne 
(2008) predicted that parsec-scale starbursts would be associ- 
ated with lower luminosity AGNs. Kawakatu & Wada (2008) 
predicted that the ratio of the scatter in the AGN to starburst 
luminosity would increase in low-luminosity AGNs. To re- 
solve this issue, it is necessary to study Seyfert galaxies over a 
wide AGN luminosity range. Thus, we performed K- and L- 
band spectroscopy of a large sample of Seyfert galaxies with 
luminosity ranging from low to high. Throughout the paper, 
H = 75km s^Mpc -1 , ft M =0.3, il A =0.7 are adopted. 

2. Target Selection 

Imanishi (2003), Imanishi & Alonso-Herrero (2004), and 
Imanishi & Wada (2004) studied Seyfert galaxies taken from 
the CfA (Huchra & Burg 1992) and 12 //m (Rush et al. 1993) 
samples. The CfA and 12 /im samples were selected through 
optical spectroscopy of large numbers of galaxies limited by 
optical and IRAS 12 /zm fluxes, respectively. The CfA sam- 
ple is usually regarded as a complete sample of optically se- 
lected Seyfert galaxies. However, Ho & Ulvestad (2001) ar- 
gued that the CfA sample is likely to be complete only for rel- 
atively bright objects. 



Therefore we added objects from the Palomar sample (Ho 
et al. 1995) in this paper. Seyfert galaxies in the Palomar 
sample are selected through optical spectroscopy of nearby, 
bright (Bt <12.5 mag), northern (declination >0°) galaxies 
covering a wider AGN luminosity range. Our aim is to con- 
firm the relationship between AGN and nuclear starburst ac- 
tivity over a wide AGN luminosity range, paying particular at- 
tention to low-luminosity AGN. To be easily observable from 
Mauna Kea, Hawaii (our observing site; latitude ^20°), the 
declinations of Seyfert galaxies were limited to greater than 
—30°. Owing to the telescope limit on the IRTF 3.0-m tele- 
scope used in this study, a restriction on declination of less 
than 68° was also applied. In total, 13 objects in the Palomar 
Sample (Table 1) were observed. We also added some addi- 
tional re-observed Seyfert galaxies to our targets. These galax- 
ies are those studied by Imanishi (2003), Imanishi & Alonso- 
Herrero (2004) and Imanishi & Wada (2004), in which PAH 
emission was barely discernible, but with less than 3cr, as re- 
ported in these papers (NGC262, NGC513, Mrk993, MCG-3- 
58-7, F03450+0055, NGC931 and NGC5548). Moreover, we 
also included 0152+06 in the CfA sample and MCG-2-8-39 in 
the 12 /im sample in our targets. Both of them meet the decli- 
nation criteria of these paper but were not observed in previous 
papers. To summarize, 13 objects from the Palomar sample, 
seven re-observed objects discussed in previous papers, and 
two new objects from CfA and 12 /im samples are included 
in this paper. 

Six of the 22 targets are Seyfert 1 galaxies and the remain- 
ing 16 targets are Seyfert 2 galaxies. Infrared (8-1000 /im) 
luminosity of our targets is logL/ii=42.5-44.8 ergs s _1 for the 
Palomar sample and logLijj=43.6-44.E ergs s _1 for the CfA 
and 12/im samples (see a caption of Table 1 for the definition 
of Lin). Detailed information on the targets is summarized in 
Table 1. We define type 1-1.5 and type 1.8-2 Seyfert galax- 
ies as Seyfert 1 and Seyfert 2 galaxies, respectively. Mrk993 
was classified as a Seyfert 2 galaxy in the CfA sample, but 
later reclassified as a Seyfert 1 .5 galaxy (Osterbrock & Martel 
1993). We adopt the latter classification and thus classified it 
as a Seyfert 1 galaxy, although it was regarded as a Seyfert 2 
galaxy by Imanishi & Alonso-Herrero (2004). 

3. Observations & Data Analysis 

Table 2 summarizes our observation log. All targets were ob- 
served with SpeX (Rayner et al. 2003) on the IRTF 3.0-m tele- 
scope on Mauna Kea, Hawaii. SpeX has a 1.9-4.2/im cross- 
dispersed spectroscopic mode (LXD1.9 mode), enabling us to 
observe K- and L-band spectra simultaneously. Because the 
seeing in K measured with the SpeX guiding/imaging mode 
was about 0."4-0."8 (FWHM) throughout the observations (ex- 
cept on 26 and 27 March 2009), we consistently used a 0"8- 
wide slit. On 26 and 27 March 2009, the seeing was over 
l'.'O at K, so we used a l'.'6-wide slit during those observa- 
tions. A standard telescope nodding technique (ABBA pattern) 
with a throw of 7'.'5 was employed along the slit to subtract 
background emission. The physical scale probed by our slit 
spectroscopy is 24-500 pc taking into account the redshift, z 
= 0.002-0.032, and slit width. We observed A-, F- and G- 
type main-sequence stars (Table 2) as standard stars. These 
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were used to correct for transmission through the Earth's at- 
mosphere. We observed the standard stars at airmasses similar 
to the targets (<0.1 difference), before and after the observa- 
tion of each target to better correct for possible time variations 
in the Earth atmosphere. 

The spectra were reduced using standard IRAF tasks. 2 
Frames taken with an A (B) beam were subtracted from frames 
taken with a B (A) beam, and then median images were taken 
and divided by a spectroscopic flat image. Known bad pix- 
els and remaining pixels hit by cosmic rays were replaced 
with interpolated values from the surrounding pixels. The 
Seyfert galaxies and standard-star spectra were then extracted 
by integrating signals over V.'3-2'.'5 along the slit. The one- 
dimensional spectrum was extracted by fitting the signals in 
each spectral order in the two-dimensional image. We used 
the wavelength-dependent nature of transmission through the 
Earth's atmosphere for our wavelength calibration. The ob- 
ject spectra were divided by spectra of standard stars to re- 
move the effects of atmospheric absorption and then multiplied 
by the spectra of blackbodies with temperatures corresponding 
to those of the individual standard stars (Table 2). Flux cali- 
bration was performed based on the signal detected inside our 
slit, taking into account the differences in exposure time be- 
tween our targets and the standard stars. K- and L-band mag- 
nitudes for the standard stars were estimated from their V-band 
magnitudes, adopting appropriate V — K and K — L colors for 
each individual standard star stellar type (Tokunaga 2000). A 
strong methane absorption line appears X res t = 3.315/mi in the 
Earth's atmosphere. Although this could be corrected by stan- 
dard stars in theory, it introduces a relatively strong amount of 
noise. Thus, we excluded a few data points at around X res t = 
3.315 pm, whose transmission were less than 40% of nearby 
signals without being affected by the absorption line. 

4. Results 

4.1. L-band spectrum 

Figure 1 shows flux-calibrated 3.0-3.6 /tin slit spectra in the 
L-band. Although the full L-band (2.8-4.2 /im) wavelength 
range is covered by SpeX, both ends of the L-band are noisy 
because of poor transmission through the Earth's atmosphere 
(A < 3.0/im) and high background emission from the Earth's 
atmosphere (A > 3.6/im). We removed those data points for 
the discussions in this paper because they were not necessary 
for our scientific aims. In these plots, some objects show a 
clear excess at 3.29 /zmx(l+z) whose wavelength corresponds 
to the redshifted 3.3 pm PAH emission feature. We regard the 
features as detected when at least two consecutive data points 
at around 3.29 /imx(l+z) are higher than 1.5 times the scat- 
ter of the continuum level. Using this criterion, objects with 
detectable PAH emission are marked with "3.3 /im PAH", and 
PAH non-detected objects are marked with "3.3 /im PAH(?)" 
in the figure. We fitted the emission line to a Gaussian func- 
tion by setting the normalization, full width of half maximum 
(FWHM) and central wavelength of the emission as free pa- 
rameters, and then estimated the flux {J^.^pah), luminosity 

2 IRAF is distributed by the National Optical Astronomy Observatories, o] 
Inc., under cooperative agreement with the N ational Science Foundation. 



(L3.3PAH) and rest-frame equivalent width (EW3.3PAH) of 
the 3.3 /im PAH emission. For objects with no clear PAH emis- 
sion features, we estimated upper limits by fitting the emission 
line with a Gaussian function whose height, central wavelength 
and FWHM were fixed as 3 a, where 1 a was a dispersion 
of its continuum emission, 3.29x(l+z)/im and 0.02x(l+z)//m 
(Tokunaga et al. 1991), respectively. The strengths of the 3.3 
/mi PAH emission features are summarized in Table 3. 

The 3.3 /im PAH emission feature is detected in ~ 23% 
(5/22) of the observed Seyfert galaxies (NGC6764, 5194, 
5273, F03450+0055 and MCG-3-58-7). NGC6764 shows a 
strong PAH emission feature, and its rest-frame equivalent 
width (EW3.3PAH) is as high as those of starburst galaxies 
(^100 nm; Imanishi & Dudley 2000). Although four objects 
(NGC5194, NGC5273, F03450+0055 and MCG-3-58-7) show 
detectable 3.3 /im PAH emission, the EW3.3PAH values are 
significantly smaller than ~400 nm. The other 17 objects do 
not display detectable PAH features. 

4.2. K-band spectrum 

Flux-calibrated 2.0-2.5 pm slit spectra in the K-band are 
shown in Figure 2. Most of the spectra show a flux depression 
at X bs > 2.3pm in the observed frame, which is attributed to 
the CO (Av=2-0) molecular absorption feature. The strength 
of the CO absorption feature was quantified by a spectroscopic 
CO spec index defined by Doyon et al. (1994) as follows: 



CO spec = -2.5log 10 (R 2 .3e), 



where (-R2.36) is an average of actual signal at A res t = 2.3 1— 
2.40 /im divided by a power-law continuum (F\ = a x A* 3 ) 
extrapolated from shorter wavelengths. As some strong emis- 
sion features appear at around 2/im in K-band, such as H2 
1-0 S(l) (A rest = 2.122/im), Br7 (A rest = 2.166/mi) and H 2 
1-0 S(0) (A rest = 2.223,um), we used data points at A rest = 
2.1 — 2.29/im to fit the continuum level after excluding these 
emission lines. The continuum level is shown as a solid line in 
Figure 2. The CO spec of all our targets are given in Table 3. 
The CO S p ec values of cool stars are typically 0.2-0.3 (Doyon 
et al. 1994). For objects with a clear CO absorption fea- 
ture (CO S p ec > 0.02), "CO absorption" signature is written 
in Figure 2. About 77 % (17/22) of observed Seyfert galax- 
ies show clear CO absorption features. Figure 3 is a histogram 
of the COspec of Luminous Infrared Galaxies (LIRGs: their 
luminosities are generally explained by starbursts) measured 
by Ridgway et al. (1994). It suggests that their typical value is 
also ~ 0.2-0.3. If the properties of nuclear starbursts in Seyfert 
galaxies are similar to those of starbursts in LIRGs, the CO spec 
of the starbursts in Seyfert galaxies is also expected to be 0.2- 
0.3. However, the observed CO spe c is smaller than the typical 
value for a LIRG because featureless AGN emission is super- 
posed in the K-band spectrum. 

ated by the Association of Univers ities for Research in Astronomy (AURA), 
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5. Discussion 

5.1. K — L color 

Because we took K- and L-band spectra simultaneously, pos- 
sible slit loss caused by IRTF/SpeX tracking error is simi- 
lar and sky conditions are the same for the K- and L-bands. 
Therefore, the derived K — L colors are very reliable. The 
emission from normal stars (>2000K) in the L-band is much 
weaker than in the K-band, whereas hot (^lOOOK) dust heated 
by an AGN emits strongly in both the K- and L-bands. Thus the 
K — L colors become smaller (bluer) with increasing relative 
contribution from stellar emission to the total emission of nu- 
clear region. Indeed, intrinsic K — L colors of Seyfert 1 galax- 
ies are typically 1-2 mag (Alonso-Herrero et al. 2003) 3 , and 
those of normal spiral galaxies are K — L < 0.4 (Willner et al. 
1984). Therefore, the color enables us to distinguish starburst- 
important Seyfert galaxies from AGN-dominant ones. The nu- 
clear K — L colors derived from our slit spectra are summarized 
in column (6) of Table 3. 

For the combined sample used in this paper, taken from 
Imanishi (2003), Imanishi & Alonso-Herrero (2004) and 
Imanishi & Wada (2004), we compared the K — L colors of 
Seyfert 1 and 2 galaxies (Figure 4). Most of the Seyfert 1 galax- 
ies are distributed in the K — L = 1-2 mag range, whereas the 
colors of almost half of Seyfert 2 galaxies are bluer than the 
Seyfert 1 galaxies (K — L < 1 mag). It shows that stellar con- 
tamination is relatively larger in Seyfert 2s than Seyfert Is in 
the observed total emission. Two possibilities could explain the 
result. First, AGN emission is strongly affected by the extinc- 
tion of dusty torus compared with nuclear starburst emission. 
In the AGN unified model, as the AGN emission comes from 
the central region, the emission from a Seyfert 2 galaxy passes 
through the dusty torus and is more strongly absorbed by the 
dusty torus than those from a Seyfert 1 galaxy. So, observed 
AGN emission of Seyfert 2 galaxies would be relatively weak 
in comparison to Seyfert 1 galaxies. Whereas, if the nuclear 
starbursts are occurring in the outer region of the dusty torus, 
then the absorption effect would not be much difference be- 
tween these types of Seyfert galaxies. In this case, the starburst 
contribution to the total emission would be relatively signif- 
icant in Seyfert 2s when compared with Seyfert Is and then 
the K — L colors of Seyfert 2s would be bluer than those of 
Seyfert Is. Another possible explanation is the starbursts in 
the central region in Seyfert 2 galaxies could be more active 
than those in Seyfert 1 galaxies. If strong nuclear starbursts in- 
flate the dusty torus more thickly than weaker starbursts, then a 
Seyfert galaxy with active starbursts is much more likely to be 
observed as type 2, and the color of it would be bluer than one 
with weak starbursts. However, as we will discuss in §6.0.1, 
we find no clear difference between Seyfert 1 and 2 galaxies in 
quantitatively measured 3.3 fim PAH luminosity. So, it is un- 
likely that the difference of the K — L colors between Seyfert 
Is and 2s is caused by the intrinsic difference of starburst ac- 
tivities in them. Therefore we suggest that the bluer K — L 
colors of Seyfert 2s than Seyfert Is are due to the obscuration 
of AGN emission and to relative contamination of the stellar 
emission in Seyfert 2 galaxies. Kotilainen et al. (1992) and 

3 K' was used in Alonso-Herrero et al. 2003; we assumed the differ ence be 



Alonso-Herrero et al. (1996) deconvolved surface brightness 
profile of Seyfert Is and 2s within 3-arcsec aperture into a com- 
bination of a nuclear source (non-stellar component), a bulge 
component and a disk component, and suggested that the stel- 
lar contamination to the total emission is larger in Seyfert 2s 
than in Seyfert Is. 

Figures 5a and 5b show the K — L color versus the 
EW3.3PAH and CO spec of the sample. These figures show 
that both the EW3.3PAH and CO spec decrease with redden- 
ing K — L color. Thus they are useful tracers for estimating 
the starburst contribution to the nuclear spectra inside our slit. 
The trends of Seyfert 2 galaxies in Figure 5 are similar to those 
for Seyfert 1 galaxies. In Figure 5b, however, some fraction 
of Seyfert 2 galaxies in our sample reach to CO spec = 0.2-0.3, 
which is closer to those of LIRGs, in comparison to Seyfert 1 
galaxies. This is good agreement with the suggestion above. In 
other words, the CO spec of some fraction of Seyfert 2 galaxies 
are not strongly diluted by AGN emission. 

6. Comparison of nuclear starbursts with AGN power 

6.0.7. IRAS 12 and 25 fim luminosity 

We compare the nuclear starburst activity with AGN activ- 
ity over a wide AGN luminosity range. The nuclear starburst 
activity is reasonably quantifiable from the observed 3.3 /im 
PAH emission luminosity inside our slit spectrum. Meanwhile 
Alonso-Herrero et al. (2003) showed that the obscuring effects 
of the AGN emission for Seyfert galaxies became insignificant 
at longer than 10 fim. IRAS 60 and 100 /im luminosities 
contain more contamination of emission from star formation 
in the host galaxy than do IRAS 12 and 25 /im luminosi- 
ties (Spinoglio & Malkan 1989; Rodriguez Espinosa & Perez 
Garcia 1997; Alonso-Herrero et al. 2002), so we use IRAS 
12 and 25 /mi luminos ities as good tracers of AGN power in 
Seyfert galaxies. A comparison between the IRAS 12 and 25 
/im luminosities and the observed nuclear 3.3 /im PAH emis- 
sion luminosity is shown in Figures 6a and 6b, respectively. 
In these figures, we combined with the data of four sources ob- 
served by Rodriguez- Ardila & Viegas (2003) using the 0."8 slit 
width of IRTF/SpeX. 

We see no clear difference in the L33PAH between the 
two types of Seyfert galaxies when normalized to the AGN 
power. This means that we see no evidence of the possibil- 
ity that Seyfert 2 galaxies tend to show intrinsically stronger 
starburst activities than do Seyfert 1 galaxies, a possibility 
discussed in §5.1. In Figure 6, the nuclear 3.3 /im PAH lu- 
minosities for our objects combined with previous data are 
L3.3PAH ~ 10 38 -10 42 ergs s" 1 . In starburst-dominated galax- 
ies, the 3.3 fim PAH-to-far-infrared (40-500 /im) luminosity 
ratios (L3,3 P ah/Lfir) are ~ 1 x 10~ 3 (Mouri et al. 1990). 
The nuclear star-formation rates (SFR) are calculated by SFR 
M Q yr" 1 = (Lfir/2.2 x 10 43 ergs s" 1 ) (Kennicutt 1998). So 
the SFR of our sample could be up to 4.5^-50 xl0~ 3 M Q yr -1 . 
Wada et al. (2009) showed that a supernova rate (SNR) of 5.4 
to 540 x 10" 5 yr _1 , corresponding to 7.7-770 x 1O~ 3 M yr -1 
(SNR ~ 0.007 SFR), can produce torus heights of as large as 
10 pc at the outer side ~ 5 pc. Hence, the star-formation rate of 

;en K' and K was negligible. 
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almost our entire sample could be high enough to create a geo- 
metrically thick dusty torus. We should note, however, that we 
have only upper limits on the nuclear 3.3 fim PAH luminosities 
for many objects in our sample, so it is likely that the sample 
contains Seyfert galaxies that do not have a thick dusty torus. 

The figures show that the nuclear starburst luminosity de- 
creases with decreasing AGN power and that the trend dose 
not change over a wide AGN luminosity range. The ratios of 
starburst activity to AGN activities for both Seyfert 1 and 2 
galaxies show the same levels of scatter and we find no obvious 
difference between them. We apply the generalized Kendall 
rank correlation statistics provided in the Astronomy Survival 
Analysis package (ASURV: Isobe et al. 1986) (which handles 
data with upper limits) to Figures 6a and 6b. The probabilities 
that a correlation is not present for the figures are both found 
to be ^0.00%. The results show that there are tight correlation 
between 3.3 PAH luminosity in central region and IRAS 
12 and 25 fim luminosities. 
6.0.2. Ground-based N-band luminosity 

Although the contamination from star formation activity in 
IRAS 12 and 25 fim is less than in the IRAS 60 and 100 fim, 
it could be that the IRAS 12 and 25 fim data still contain a sig- 
nificant amount of emission from star-forming activity of host 
galaxies in our low luminosity AGN sample, given the large 
aperture of IRAS 12 fim (075 x 4.'5) and 25 fim (0.75 x 4f6). 
Ramos Almeida et al. (2009) compared small aperture unre- 
solved nuclear 10 fim emission (0'.'4-0'.'5) and 20 fim emis- 
sion (0''5-0"6) with the IRAS 12 fim and 25 fim emission, 
respectively, and found the large aperture data are largely con- 
taminated by starlight. So it is likely that the IRAS 12 fim and 
25 fim emission is significantly contaminated by star formation 
in host galaxies. To reduce the possible ambiguity, we should 
compare the 3.3 fim PAH emission with an AGN indicator 
measured with a smaller aperture. The N-band (Ao=10.78/im, 
AA=5.7 fim) is observable from the ground, and its luminosity 
is also thought to be a good tracer of AGN power (Alonso- 
Herrero et al. 2002; Ramos Almeida et al. 2009; Levenson et 
al. 2009). Gorjian et al. (2004) presented N-band photometric 
data for the central regions of Seyfert galaxies. They used the 
Palomar 5-m telescope with a 1"5 aperture. Only 41% (9/22) 
of our objects were measured by Gorjian et al. (2004). We 
combine these with data from Imanishi & Wada (2004), which 
is shown in Figure 7. Although the number of low-luminosity 
AGNs is small and their 3.3 PAH emission luminosities 
are only upper limits, the three low luminosity Seyfert galax- 
ies are distributed around a line extrapolated from the stronger 
AGN activity region. We also apply the generalized Kendall 
rank correlation statistics to Figure 7. The uncorrected prob- 
ability is ^0.03% for the figure. These results mean that the 
correlation between the luminosities of the nuclear starbursts 
detected inside the slit spectra and central AGNs is statistically 
confirmed in Seyfert galaxies for all observed quantities. 

We found the relation between the luminosities of 3.3 fim 
PAH emission within our slit width and of N-band. However 
it is worried that the L^^pah is affected by physical aperture 
size probed by the slits. It means that we are concerned about 
the possibility that the L^^pah of the most luminous object is 
due to the largest area we observed. If star formation traced by 
L3.3PAH is spatially extending in a wide area of host galaxy, 



then the L^^pah becomes higher with lager physical scale. 
To explore this, we compare the L3 .spah with physical area 
of each source (Figure 8). The figure clearly shows that not all 
high L3.3PAH objects are wide physical area coverage. That 
is to say, the star formation would be occurring in the central 
region inside the slits and not spatially extend. 

6.1. Trigger of accretion onto SMBH 

We find no evidence that the ratio of starbursts to AGN lumi- 
nosity deviates upward in low-luminosity AGNs, as predicted 
by Kawakatu & Wada (2008), or that nuclear starburst activ- 
ities increase with decreasing AGN activities as predicted by 
Ballantyne (2008). Furthermore, we find no clear change in 
the correlations over a wide AGN luminosity range. This is ex- 
plained if the main mechanism that connects the AGN and nu- 
clear starbursts is the same at each luminosity range. A model 
shown in Wada & Norman (2002; see also Wada et al. 2009), 
which predicts the enhancement of a mass accretion rate onto 
a central SMBH owing to increased turbulence of molecular 
gas in the torus caused by nuclear starbursts, is one possible 
scenario that could account for the luminosity correlation. 

6.2. The origin of the CO spec index 

We calculated stellar luminosity (Lk stellar) from CO spec . 
If nuclear starbursts with properties similar to star-formation- 
dominated LIRGs appear inside our slit and the entire K-band 
emission comes from stars, then the detected CO spec must be 
0.2-0.3. We assumed the original CO spec of the cool star was 
0.25. When the contributions of stellar emission to total emis- 
sion are 100%, the average signal at A res t = 2.31-2 .4 fim is re- 
duced by about 20.5% compared with an extrapolation from a 
shorter wavelength. When the contribution of AGN emission to 
the total emission is 100 %, the absorption feature does not ap- 
pear. When the contributions of both stellar and AGN emission 
are equal, the average signal at A rest = 2 31-2 A fim should be 
10.25% reduced from the shorter wavelength, and the CO spec 
should be 0.1 17. In the case that the CO spec is 0.1 or 0.15, the 
contributions of stellar emission to total emission should be 
43% or 63%, respectively. So, we can estimate the Lxsteiiar- 
Figure 9 compares the Lk stellar and the 3.3 fim PAH emission 
luminosity in the nuclear region detected inside our slit spectra 
combined with Fig. 5 of Imanishi & Alonso-Herrero (2004). 
Although they plotted these results for Seyfert 2 galaxies only, 
we included in Figure 9 Seyfert 1 galaxies whose CO spe c were 
derived in Imanishi & Wada (2004). 

The K-band to infrared luminosity ratios of starburst- 
dominated LIRGs were estimated to be Lk / £7^(8-1000^111) ~ 
10 -i.6±o.2 (Go i dader et a i 1997) If the nuc i ear starbursts 
in Seyfert galaxies have properties similar to the starburst- 
dominated LIRGs, then the same relation should hold. The 
3.3 fim PAH to infrared luminosity ratios in starbursts were 
found to be L 3 . SPAH / '£f/fl(40-500/«ii) ~ 10~ 3 (Mouri et al. 
1990). Assuming Lip ~ Lpip for starbursts, the luminosity 
of the 3.3 fim PAH emission to the K-band stellar luminosity 
ratios are expected to be L^^pah / LKsteiiar ~ 10 -14 , if both 
luminosities trace the same starbursts. The value is shown as 
a solid line in Figure 9. Most of our Seyfert 1 galaxies and 
about half of the Seyfert 2 galaxies plotted in the figure are dis- 
tributed around the line. However, the remaining half of the 
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Seyfert 2 galaxies are distributed under the line. Given that 
in Figure 6, the L^.?,pah distribution of both Seyfert 1 and 2 
galaxies is similar within their dispersions and no obvious dif- 
ference is seen between them, the discrepancy must be caused 
by L K stellar- As we discussed above, the L KsteUar is esti- 
mated from the CO spec value, which is related to the fraction 
of stellar to AGN emission. Nuclear starbursts are expected 
to occur at the outer part of the dusty torus and the AGN ra- 
diation of type 2 Seyfert galaxy is significantly attenuated by 
the dusty torus. Therefore, the CO sp ec is likely to increase 
and thus increase LKsteUar in Seyfert 2 galaxies. In contrast, 
because L^^pah is estimated from the flux of 3.3 /im PAH 
emission, PAH-estimated starburst luminosity is not affected 
by the absorption of AGN emission. In summary, we suggest 
that because LKsteUar represents a relative value of stellar to 
AGN emission and L^^pah represents an absolute value of 
stellar emission, it is possible that LKsteUar is overluminous 
and that the ratio of Ls.spah / Lxstellar decreases in Seyfert 2 
galaxies. 

The other possibility is that the signatures of stellar emission 
detected in the K-band spectra are significantly contaminated 
by old bulge stars in the nuclear part of the host galaxy, whereas 
the 3.3 /im PAH emission originates in spatially unresolved 
(smaller than subarcsecond) nuclear starbursts in the dusty 
torus because old stars do not have enough PAH-exciting UV 
photons. Ivanov et al. (2000) and Imanishi & Alonso-Herrero 
(2004) suggested that the emission of old stars in the host 
galaxy mainly produces the K-band spectra. However, most 
of the Seyfert 1 and some fraction of the Seyfert 2 galaxies 
have L3.3P ah I 'LKsteUar ratios similar to the value expected 
when the origin of both luminosities is the same. Therefore, 
although the old stars in host galaxy may contribute to the ob- 
served K-band spectra of some fraction of Seyfert galaxies, it 
is not likely that old stars dominate the K-band spectra in all 
objects. 

7. Conclusion 

The results of an infrared K- and L-band spectroscopic study 
of 13 Seyfert galaxies from the Palomar sample, and seven 
Seyfert galaxies with non-detectable 3.3 /im PAH emission dis- 
cussed in previous papers and two new Seyfert galaxies from 
CfA and 12/im samples are presented. Our L-band spectro- 
scopic method successfully detected the 3.3 /im PAH emission 
in ~ 23 % (=5/22) of the observed Seyfert galaxies. Also, 
our K-band spectroscopic method showed ~ 77 % (=17/22) of 
Seyfert nuclei in our sample have clear CO absorption feature. 
We examined the relationship between AGN activity and nu- 
clear starburst activity over a wide AGN luminosity range using 
our spectra together on previously published spectra (Imanishi 
2003; Imanishi & Alonso-Herrero 2004; Imanishi & Wada 
2004). Our conclusions are summarized as follows: 

1. The K — L colors of Seyfert 2 galaxies are widely dis- 
tributed toward the blue, compared with those of Seyfert 
1 galaxies. This implies either that the dusty torus ab- 
sorbs the AGN emission of Seyfert 2 galaxies from the 
central region or that Seyfert 2 galaxies tend to have 
stronger nuclear starbursts than do Seyfert 1 galaxies. 



As the L3 3PAH of Seyfert 1 and 2 galaxies do not differ 
significantly, we have interpreted this as the result of the 
effect that the AGN emission of Seyfert 2s is absorbed 
by dust of torus and then the stellar emission becomes 
relatively larger. 

2. The L3.3PAH shows the same range in Seyfert 1 and 2 
galaxies. The star-formation rates of our sample are up 
to 4.5-450 xl0~ 3 M Q yr -1 , which could be enough to 
swell the dusty torus by supernova explosion; in short, 
it is geometrically thick, although many objects we pro- 
vided are only with upper limits on the 3.3 fim PAH lu- 
minosities. 

3. The L3.3PAH correlates with mid-infrared (N-band) lu- 
minosity with small aperture (= tracing AGN activities), 
and their luminosity ratio does not vary significantly over 
a wide AGN luminosity range. Moreover the L^^pah 
is independent of physical scale we probed with slits. 
Therefore the L3 .zpah (i-e. star formation) would con- 
centrate on the central region, and the nuclear star for- 
mation would induce the accretion onto the SMBH and 
encourage the AGN activity. In this work, we find no 
evidence that nuclear starbursts are stronger in lower- 
luminosity AGN. 

4. We suggest that the 3.3 /im PAH emission luminos- 
ity and K-band stellar luminosity (LKsteUar) originate 
in the same phenomenon. However, for some fraction 
of Seyfert 2 galaxies, LKsteUar is likely to be overes- 
timated because of flux attenuation of AGN emission 
caused by a dusty torus and contaminated by old stars 
in the spheroid of the host galaxy. 
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Table 1. Target information 



Object 


z 


D L 


Class 


fl2 


/25 


.ho 


/ioo 


/l0.8 


logLi R 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Palomar sample 


NGC676 


0.0050 


20.2 


Sy2 


0.04 


0.06 


0.27 


0.80 




42.3 


NGC1167 


0.0165 


66.8 


Sy2 


<0.03 


<0.04 


0.12 


1.05 




43.3 


NGC1275 


0.0176 


71.2 


Syl 


0.93 


3.02 


7.09 


7.60 




44.8 


NGC1358 


0.0134 


54.3 


Sy2 


<0.08 


<0.12 


0.38 


0.93 




43.4 


NGC3982 b ' c 


0.0035 


14.2 


Sy2 


0.47 


0.97 


7.18 


16.20 


<0.1 1 


43.3 


NGC4258 6 


0.0015 


6.2 


Sy2 


2.25 


2.81 


2.81 


78.40 




43.0 


NGC4579 b ' c 


0.0050 


20.3 


Sy2 


1.12 


0.78 


5.93 


21.40 


0.08(0. 06*) 


43.7 


NGC5194 & ' C 


0.0020 


8.0 


Sy2 


7.21 


9.56 


97.40 


221.00 


<0.0 7 


43.9 


NGC5273 b 


0.0036 


14.4 


Syl 


0.12 


0.29 


0.90 


1.560 




42.5 


NGC6764 


0.0080 


32.4 


Sy2 


0.54 


1.33 


6.62 


12.40 




44.0 


NGC6951 


0.0048 


19.1 


Sy2 


1.34 


2.16 


16.20 


41.80 




44.0 


NGC7479 6 


0.0079 


32.0 


Sy2 


1.37 


3.86 


14.90 


26.70 




AAA 


NGC7743 b 


0.0057 


22.9 


Sy2 


0.10 


0.18 


0.92 


3.40 




43.0 


CfA & 12 ^m 


NGC262 C 


0.0151 


61.3 


Sy2 


0.31 


0.84 


1.29 


1.55 


0.12 


44.0 


NGC513 C 


0.0195 


79.4 


Sy2 


0.17 


0.28 


1.94 


4.05 


<0.10 


44.2 


NGC93F 


0.0167 


67.5 


Syl 


0.61 


0.61 


1.32 


4.55 


0.21 


44.3 


NGC5548 Q ^ C 


0.0172 


67.6 


Syl 


0.40 


0.77 


1.07 


1.61 


0.29 


44.1 


Mrk993 fc ' c 


0.0155 


62.8 


Syl 


<0.13 


<0.13 


0.30 


1.32 




43.6 


MCG-2-8-39 c 


0.0299 


122.3 


Sy2 


0.20 


0.48 


0.51 


<1.42 




AAA 


MCG-3-58-7 c 


0.0320 


131.1 


Sy2 


0.28 


0.80 


2.42 


3.36 


0.20 


44.8 


F03450+0055 c 


0.0310 


126.9 


Syl 


0.28 


0.51 


0.47 


<3.24 


0.10 


44.6 


0152+06 b 


0.0174 


70.6 


Sy2 


<0.08 


<0.16 


0.50 


1.15 




43.7 



Col.(l): object name, a, b, c are reference; a)Palomar sample (Ho & Ulvestad 2001), b)CfA (Huchra & Burg 1992), c)12 
/j,m (Rush et al. 1993). Col.(2): redshift. Col.(3): luminosity distance in Mpc. Col.(4): AGN class. Col.(5)-(8): flux at 12/im, 
25/im, 60/im, and 100^m in Jy, respectively from IRAS faint source catalog. Col. (9): ground-based 10.8/im photometric data, 
observed with the Palomar 5-m telescope using a l."5 aperture (Gorjian et al. 2004). NGC4579 was observed also by Horst et al. 
(2009) using VLT/VISIR with a l."27 small aperture. Col. (10): decimal logarithm of infrared (8-1000/iin) luminosity in ergs s _1 
calculated with L IR = 2.17 x 10 39 x L> L (Mpc) 2 x(13.48x/i 2 +5.16x/ 2 5+2.58x/ 60 +/ioo) ergs s" 1 (Sanders & Mirabel 1996). 



No. ] Oi, Imanishi & Imase 9 



Table 2. Observing Log 



Object 


Date 


Itime 


Slit 


Std star 


Type 


V 


K 


L 


T eff 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Palomar sample 


NGC676 


2009 Oct. 27 


48 


0.8 


HD 13043 


G2V 


6.87 


5.39 


5.34 


5830 


NGC1167 


2009 Oct. 26 


70 


1.6 


HD 19600 


A0V 


6.42 


6.32 


6.32 


9480 


NGC1275 


2009 Oct. 27 


44 


0.8 


HD20995 


A0V 


5.61 


5.61 


5.61 


9480 


NGC1358 


2009 Oct. 28 


120 


0.8 


HD21019 


G2V 


6.20 


4.42 


4.37 


5830 


NGC3982 


2005 Apr. 29 


68 


0.8 


HR4496 


G8V 


5.33 


3.53 


3.47 


5430 


NGC4258 


2005 Apr. 27 


80 


0.8 


HR4572 


F6V 


6.62 


5.41 


5.37 


6385 


NGC4579 


2005 Apr. 28 


40 


0.8 


HR4708 


F8V 


6.40 


5.05 


5.01 


6130 


NGC5194 


2005 May. 1 


60 


0.8 


HR4845 


GOV 


5.95 


4.54 


4.49 


5930 


NGC5273 


2005 Apr. 29 


60 


0.8 


HR4845 


GOV 


5.95 


4.54 


4.49 


5930 


NGC6764 


2005 May. 1 


60 


0.8 


HR7294 


G4V 


6.57 


5.04 


4.99 


5740 


NGC6951 


2005 Apr. 30 


100 


0.8 


HR7783 


G3V 


5.93 


4.44 


4.39 


5785 


NGC7479 


2009 Oct. 26 


56 


0.8 


HD217577 


G2V 


8.66 


7.12 


7.07 


5830 


NGC7743 


2009 Oct. 27 


75 


1.6 


HD377 


G2V 


7.60 


6.12 


6.07 


5830 


CfA & 12 fxm 


NGC262 


2007 Aug. 31 


80 


0.8 


HR410 


F7V 


6.31 


4.99 


4.95 


6240 


NGC513 


2007 Aug. 27 


120 


0.8 


HR410 


F7V 


6.31 


4.99 


4.95 


6240 


NGC931 


2007 Aug. 26 


96 


0.8 


HR720 


GOV 


5.89 


4.48 


4.43 


5930 


NGC5548 


2005 Apr. 29 


20 


0.8 


HR5346 


F8V 


6.25 


4.90 


4.85 


6130 


Mrk993 


2007 Aug. 31 


112 


0.8 


HR410 


F7V 


6.31 


4.99 


4.95 


6240 


MCG-2-8-39 


2007 Aug. 30 


136 


0.8 


HR784 


F6V 


5.78 


4.57 


4.53 


6385 


MCG-3-58-7 


2007 Aug. 30 


80 


0.8 


HR8457 


F6V 


6.09 


4.88 


4.84 


6385 


F03450+0055 


2007 Aug. 28 


112 


0.8 


HR962 


F8V 


5.06 


3.71 


3.67 


6130 










HR145 


F7V 


6.41 


5.09 


5.05 


6240 


0152+06 


2007 Aug. 29 


120 


0.8 


HR508 


G3Va 


6.27 


4.78 


4.73 


5785 



Col.(l): object name. Col. (2): observing date in UT. Col. (3): net on-source exposure time in minutes. Col. (4): slit width in 
arcsec. Col. (5): name of standard stars. HR962 was observed as the standard star for F03450+0055. However, because it was 
saturated in some part of the K^band spectrum, we used HR145, which was observed on the same day, as the K4oand standard 
star for F03450+0055. Col. (6): stellar spectral type. Col.(7)-(9): V, K and L^band magnitude. Col.(10): effective temperature in 
K. 
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Table 3. Properties of the nuclear 3.3 PAH emission and CO absorption feature s 



Object 


/3.3PAH 


log(L 3 , 3P AH) EW 3 . 3 PAH 


CO spec 


K-L 


SFR 




(xl0 14 ergs s _1 cm -2 ) 


(ergss x ) 


(nm) 




(mag) 


(M Q yr- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Palomar sample 


NGC676 


< 1 1.58 


<39.73 


<10.79 


0.20 


-0.08 


<0.24 


NGC1167 


< 7.85 


<40.60 


< 18.02 


0.23 


0.14 


<1.80 


NGC1275 


<31.58 


<41.26 


<4.26 





1.79 


<8.22 


NGC1358 


< 9.21 


<40.49 


<21.53 


0.18 


-0.04 


<1.39 


NGC3982 


< 6.17 


<39.15 


<28.42 


0.16 


0.27 


<0.06 


NGC4258 


<13.57 


<38.77 


<4.17 


0.16 


0.95 


<0.03 


NGC4579 


< 3.22 


<39.17 


< 1.37 
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Fig. 6. (a): IRAS 12 fim luminosity, defined as j/L„ (abscissa) versus the 3.3 fim PAH emission luminosity detected inside our slit spec- 
tra (ordinate). Indices are the same as in Fig. 5. Seyfert 1 galaxies studied by other groups (IH1934-063, Mrk766, NGC3227, and 
NGC4051; Rodriguez-Ardila & Viegas 2003) are plotted as filled triangles, (b): Same as (a), but the abscissa is IRAS 25 fj,m luminosity. 
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Fig. 7. Same as Fig. 6, but the abscissa is ground-based N-band luminosity measured with small aperture. Almost all N-band data in the abscissa 
are 10.8 fim luminosity measured with a 1'.'5 aperture (Gorjian et al. 2004). Galliano et al. (2005) observed 10.4 /im luminosity with 2-arcsec 
aperture for NGC2992 and NGC7469 using ESO 3.6-m telescope/TIMMI2, and Horst et al. (2009) observed 10.49 fim luminosity with 1.27-arcsec 
aperture for NGC4579, Mrk509 and MCG -3-34-64 using VLT/VISIR. We adopted the new values for the five objects. Among the four sources 
studied by Rodriguez-Ardila & Viegas (2003), only two objects (NGC3227 and NGC4051) are plotted because of the availability of measured 
10.8 iim luminosity (Gorjian et al. 2004). The abscissa is a good indicator of AGN power, and the ordinate probes nuclear starburst luminosities. 
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Fig. 8. Relationship between physical scale we probed with slits and the 3.3 fim PAH emission luminosity inside the slits 
(same symbols as in Figure 6). Since all physical areas of objects except 5 objects (Mrk34, Mrk78, Mrk273, Mrk463, 
Mrk477) plotted in the figure are smaller than 2 square kpc, larger figure focus on the physical area from 0.001 to 1.827 
kpc 2 to show the distribution. A figure including these 5 objects are pasted small on the lower right of larger figure. 
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Fig. 9. Comparison between the K-band stellar luminosity (abscissa) and the nuclear 3.3 fim PAH luminosity detected inside our slit spectra (ordi- 
nate). Symbols are the same as in Figure 5. The length of the bar at the lower right of the figure represents the range in uncertainty of LxatMar 
due to COspec = 0.2 — 0.3. The solid line shows the predicted ratio of L3 ,3pah I ' Lpcstellar (~ 10 — 1 ) for starbursts as seen in LIRGs. 



